arXiv:1503.01735v2 [cond-mat.mes-hall] 24Jun2015 


Spin and charge transport in graphene-based spin transport devices with Co/MgO 

spin injection and spin detection electrodes 


F. Volmer, M. Drogeler, G. Giintherodt, C. Stampfer, B. Beschoten 
2nd Institute of Physics and JARA-FIT, RWTH Aachen University, D-52074 Aachen, Germany 


Abstract 

In this review we discuss spin and charge transport properties in graphene-based single-layer and few-layer spin-valve 
devices. We give an overview of challenges and recent advances in the field of device fabrication and discuss two of our 
fabrication methods in more detail which result in distinctly different device performances. In the first class of devices, 
Co/MgO electrodes are directly deposited onto graphene which results in rough MgO-to-Co interfaces and favor the 
formation of conducting pinholes throughout the MgO layer. We show that the contact resistance area product ( R C A ) 
is a benchmark for spin transport properties as it scales with the measured spin lifetime in these devices indicating 
that contact-induced spin dephasing is the bottleneck for spin transport even in devices with large R C A values. In a 
second class of devices, Co/MgO electrodes are first patterned onto a silicon substrate. Subsequently, a graphene-hBN 
heterostructure is directly transferred onto these prepatterned electrodes which provides improved interface properties. 
This is seen by a strong enhancement of both charge and spin transport properties yielding charge carrier mobilities 
exceeding 20000cm 2 /(Vs) and spin lifetimes up to 3.7 ns at room temperature. We discuss several shortcomings in 
the determination of both quantities which complicates the analysis of both extrinsic and intrinsic spin scattering 
mechanisms. Furthermore, we show that contacts can be the origin of a second charge neutrality point in gate dependent 
resistance measurements which is influenced by the quantum capacitance of the underlying graphene layer. 
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1. Introduction 

In spin-based electronics, three main aspects have to 
be considered when exploring suitable spin transport ma¬ 
terials and material combinations: (1) electrical injection 
and detection of spins, (2) their manipulation, and (3) 
the transport of spins in the material w- What makes 
graphene a promising material in the field of spintron- 
ics, is its unique spin transport performance in particu¬ 
lar at room temperature |4| where spin lifetimes of up 
to 3.7 ns [5] and spin diffusion length of 12 /xm [B] have 
been measured by means of electrical Hanle spin preces¬ 
sion measurements in non-local spin-valve devices. The 
corresponding charge carrier mobilities in these devices are 
above 20000 cm 2 /(Vs) [5j|6j. Other interesting materials 
in the field of spintronics, e.g. Si, also exhibit nanosecond 
spin lifetimes at room temperature but fail short to gra¬ 
phene in respect to the spin diffusion lengths and charge 
carrier mobilities [7]. We note that we only compare de¬ 
vice and material properties from electrical spin preces¬ 
sion measurements at room temperature. Less invasive 
spin sensitive methods such as electron spin resonance or 
optical pump-probe methods can yield much longer spin 
lifetimes, especially at low temperatures HEHg. 
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Graphene has an extraordinary band structure m and 
a weak intrinsic spin-orbit coupling at energies close to the 
Dirac point m- Often, this property is mentioned as a 
favorable aspect for graphene spintronics, because for the 
most prominent spin relaxation mechanisms the spin re¬ 
laxation rate scales with the spin-orbit coupling strength 
W- Accordingly, initial calculations promised quite long 
spin lifetimes in pristine graphene flakes up to the ms 
regime m- But the experimental values are orders of 
magnitude smaller than these predictions and only exhibit 
spin lifetimes in the range of 20 ps to 3.7 ns at room tem¬ 
perature 15, Sj IT3M39| . More elaborated calculations which 
included novel spin scattering mechanisms such as reso¬ 
nant scattering by magnetic impurities m or entangle¬ 
ment between spin and pseudospin by random spin orbit 
coupling m can explain these short spin lifetimes. Nev¬ 
ertheless, final answers about both the limiting spin relax¬ 
ation mechanism and the maximal achievable spin lifetime 
in graphene are still missing. Considering the latter, it is 
interesting to mention that electron spin resonance (ESR) 
experiments in synthesized graphene flakes yield spin life¬ 
times of conduction electrons in the range of 30 — 65 ns, de¬ 
spite of a significant defect density in the studied samples 
[[42] .43] . It has been suggested that the much longer spin 
lifetimes in ESR is due to the fact that the graphene sheets 
are free from substrate effects and metallic electrodes [113 . 
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There is a large number of publication from different 
groups about electrical injection and detection of spins in 
graphene (e.g. [51 [HI IT5H55] j. But although electric fields 
from back and/or top gates or biases along the graphene 
channel can strongly modify spin and charge transport 
properties, the actual spin precession in these experiments 
is always triggered by external magnetic fields. The rea¬ 
son for this is that the before-mentioned weak spin-orbit 
coupling in graphene is also a mixed blessing because for 
spin manipulation, e.g. via the Bychkov-Rashba effect, a 
strong spin-orbit coupling is needed which, on the 

other hand favors spin dephasing and spin scattering as it 
has been explored in III-V semiconductors [551 ; 4Sj ■ There¬ 
fore, one emerging topic in graphene-based spintronic re¬ 
search is the partial functionalization of graphene with the 
goal to achieve other ways for spin manipulation, e.g. by 
electrostatic gating. The aim is to use high quality gra¬ 
phene parts as leads for efficient spin transport whereas 
the spin manipulation is realized in a functionalized part 
of the graphene device with enhanced spin-orbit coupling. 
Routes to increase the spin-orbit coupling include the spin- 
orbit proximity effect, in which additional materials such 
as the two dimensional transition metal dichalcogenides 
(e.g. tungsten disulfide 03 ) or ferromagnetic insulators 
(e.g. EuO ^S]) are put in direct contact to graphene. 
Also hydrogenation of graphene [49] or the deposition of 
heavy adatoms [50] can yield enhanced spin-orbit coupling 
in graphene, but as a drawback, these approaches are al¬ 
ready known to strongly change its band structure [501 M- 

Very recently, an excellent review article about gra¬ 
phene spintronics compiled a significant part of both theo¬ 
retical and experimental work which has been carried out 
on spin phenomena in graphene [4]. In particular, this 
review focuses on how to measure spin transport in gra¬ 
phene, the spin-orbit coupling in pristine and modified gra¬ 
phene, magnetic moments from defects and adatoms, and 
the open question about which spin relaxation mechanisms 
limit spin transport in graphene. In our article, we expand 
the review to open questions and challenges in the experi¬ 
mental field of work. In section [2] we discuss new routes in 
device fabrication by the usage of novel transfer techniques 
and outline the issue of device contaminations during fab¬ 
rication. Then, we cover the influence of metallic con¬ 
tacts on the underlying graphene in section [3] and present 
new results on the appearance of a contact-induced sec¬ 
ond charge neutrality point in gate dependent resistance 
measurements. Next, the reliability of extracted values 
for both spin lifetimes and charge carrier mobilities is dis¬ 
cussed in sections [4] and [5] respectively. In section [5] we 
briefly comment on some experimental studies which were 
used to investigate the relevant spin relaxation mechanism 
in graphene. 

2. Methods of device fabrication 

In recent years, it has been demonstrated that elec¬ 
trical injection and detection of spins in graphene can be 


accomplished by a variety of different electrode materials, 
such as Al 2 0 3 /Co [T3HT51 [T51 ISO] . MgO/Co [51 H51 Ell 1241 - 
127] , submonolayer TiC^/MgO/Co [20] |22] [23], Ti02/Co 
d [T7j, Cu/NiFe [52], amorphous carbon/Co [53], fluori- 
nated graphene/NiFe [54], hydrogenated graphene/NiFe 
[35] , PTCA/ALD Al 2 03/NiFe [35], only Co [251123113T], 
h-BN/NiFe [32], h-BN/Co [33] [35] [35] , and YO/Co [37] , 

However, as mentioned in the introductory part, the mea¬ 
sured spin lifetimes are only in the range of 20 ps to 3.7 ns 
at room temperature. This is well below the lifetimes of 
30—65 ns measured by ESR experiments in graphene flakes 
without any contacts [52J [43]. Next to these ESR exper¬ 
iments there are also other studies which indicates that 
an insufficient barrier quality can be the bottleneck for 
the overall spin transport [201 EH [271 [55]. In the follow¬ 
ing, we therefore address the challenges and summarize 
the progress in the held of device fabrication. 

First, we discuss several shortcomings when growing 
electrode material directly onto the graphene surface. Be¬ 
cause often the starting point of a graphene-based spin 
transport device is a randomly exfoliated graphene flake 
which is typically deposited onto Si ++ /Si0 2 - The next 
fabrication step is usually a lithography process on top of 
graphene which is needed for the deposition of ferromag¬ 
netic electrodes. But as soon as graphene gets into contact 
with an organic resist, it is extremely difficult to com¬ 
pletely remove the contamination of hydrocarbons. This 
is illustrated in the atomic force microscopy (AFM) im¬ 
ages of figures [IJa)-(c), which show one of our earlier but 
not optimized dosage tests. As a resist, we used PMMA 
(950K) which was dissolved in ethyl lactate and n-butyl 
acetate with a thickness of 250 nm after spin-coating and 
baking. For developing we used isopropyl alcohol and 
methyl isobutyl ketone with a developing time of 105 s. 
E-beam writing was performed with an acceleration volt¬ 
age of 10 kV. The dosages in figures [l]( a)-(c) are given in 
fractions of 100 /iC/cm 2 . The scale in figured! a) was cho¬ 
sen to depicted the flake in high contrast. The resist is 
underdeveloped for a dosage below 70/iC/cm 2 . But also 
for dosages larger than 110 fiC/cm 2 an increasing contami¬ 
nation of the graphene flake can be observed, which can be 
easier be visualized in the phase signal of the AFM image 
in tipping mode (figure [ljb) ). The difference in contrast 
can be understood by different attenuations of the AFM 
cantilever for Si0 2 , graphene on Si0 2 , and PMMA on gra¬ 
phene. The increasing contamination for higher dosages 
can be explained by crosslinking of PMMA [55] or e-beam- 
induced defects in the graphene flake [571158] which enables 
a stronger binding with hydrocarbons. We observe the 
cleanest graphene surface at a dosage of 90/LtC/cm 2 and 
were able to reduce the amount of leftover resist residues 
even further by increasing the developing time to 210 s. 

But even for optimized lithography conditions the gra¬ 
phene flake will be contaminated with hydrocarbons on 
an atomic scale. Such contaminations are clearly seen by 
transmission electron microscopy and annealing tempera¬ 
tures of up to 700 °C are needed to remove the hydrocar- 
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Figure 1: (Color online) (a) AFM image of a graphene flake with spin coated PMMA resist. The fields were written by e-beam lithography 
with different dosages as given by the denoted factors times the reference dosage of 100/iC/cm 2 . The scale is adjusted to maximize the 
contrast of the graphene flake, (b) In the phase signal of the AFM the resist residues on top of graphene are clearly visible, (c) Image from 
(a), but now the scale is optimizes to visualize the underdeveloped resist residues, (d) AFM image of a 3 nm thick MgO layer on top of 
graphene demonstrating the Volmer-Weber island growth: rms roughness of 0.4 nm and peak-to-peak values of 2 nm. (e) Transfer technique 
for the fabrication of graphene spin-valves from the bottom-up approach in ref. [5]- Figures reproduced with permission from: (d), ref. El. 
©2014 American Physical Society; (e), ref. [5], ©2014 American Chemical Society. 


bons [55, 60 . But such high temperatures are far beyond 
the glass transition temperature of the resist and there¬ 
fore will destroy the patterning mask. Next to thermal 
annealing also other methods are applied to remove resist 
residues after the development step of the resist, e.g. low- 
density inductively coupled Ar plasma [61], CO 2 cluster 
cleaning [62], or ultraviolet-ozone treatment [63, :64j . But 
it still has to be seen if these cleaning methods are able to 
remove even the last atomic layer of hydrocarbons from the 
graphene flake. While in most of these publications only 
rather macroscopic analysis methods such as Raman spec¬ 
troscopy or even qualitative electrical measurements are 
used, only the high-resolution X-ray photoelectron spec¬ 
troscopy in case of the ultraviolet-ozone treatment of ref¬ 
erence [BH may really demonstrate the complete removal 
of all contaminations. There is, however, also a drawback 
to this method as long treatments can also create defects 
in the graphene flake. 

The next issue arises because of the chemically inert 
nature of graphene and its sp 2 hybridization, which leads 
to unfavorable growth conditions of many materials on top 
of graphene. This can already be seen in the stronger accu¬ 
mulation of the PMMA on the graphene flake compared to 
the area of Si0 2 for dosages lower than 70 /xC/cm 2 in figure 


0b) (in this figure the scale is optimized to the height of 
the underdeveloped PMMA). The reason for this accumu¬ 
lation is the low wettability of graphene [B3| and again it is 
transmission electron microscopy which reveals the cluster¬ 
ing of almost every metal deposited on graphene HHHl]. 
The interaction between metals and graphene is so weak, 
that in sub-monolayer growth of metal layers on graphene, 
the metal atoms solely resides on the hydrocarbon contam¬ 
ination and not on the clean graphene parts [601166[ . For 
example, the direct growth of MgO on graphene is gov¬ 
erned by the Volmer-Weber island growth mechanism [67] 
resulting in inhomogeneous oxide layers with presumably 
pinholes m- This is illustrated in the AFM image of a 
3 nm thick MgO layer grown on top of graphene in figure 
0d), which exhibits peak-to-peak values of up to 2 nm. 
The use of a Ti/Ti0 2 wetting layer yields more homoge¬ 
neous MgO layers Sz| but the impact of Ti on the prop¬ 
erties of graphene is still an open question [6SJ [65]. For 
Al 2 03 the evaporation of A1 and its subsequent oxidation 
to Al 2 03 can also yield rough layers with pinholes COJ- 
On the other hand, there is a report on homogeneous and 
pinhole free Al 2 03 oxide barriers on graphene by argon 
sputter deposition m- But the same group also demon¬ 
strated by Raman spectroscopy that this technique may 
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Figure 2: (Color online) (a) and (b) Gate dependent graphene resis¬ 
tance for spin transport devices at room temperature with spin life¬ 
times of 80 ps and 2 ns (data taken from reference EZD- Devices with 
long spin lifetimes exhibit a 2nd charge neutrality point at negative 
gate voltages (see red arrow in (b)). In addition to a doping of the 
underlying graphene layer, the interaction of the Co/MgO electrodes 
with graphene can suppress the electric field effect by the gate voltage 
Vq which can result in (c) pinning (devices with R C A < lkfi/im 2 ) 
or (d) no pinning (devices with R C A > lkQ/mi 2 ) of the Fermi level 
in graphene parts underneath the electrodes (corresponding Dirac 
cones are indicated by dashed lines). 


also induces defects in graphene m ■ By now, there are 
only few studies reporting epitaxial-like growth of materi¬ 
als on graphene, e.g. EuO [48j or Ti [55] . 

The hydrocarbon contaminations from the lithography 
step and the clustering of many materials on top of gra¬ 
phene are important issues for graphene-based spin trans¬ 
port devices because both can reduce the quality of the 
insulating oxide barrier between graphene and the ferro¬ 
magnetic electrodes. But as mentioned above, an insuffi¬ 
cient barrier quality can be the bottleneck for the overall 
spin transport [201 [261 HU Eg. A possible way to circum¬ 
vent some of these problems was paved by the introduction 
of mechanical transfer techniques for 2d materials (see e.g. 
E2H321). These transfer techniques provide a controlled 
deposition of flakes with high spatial precision, which re¬ 
sults in two important advantages: (1) The device fabrica¬ 
tion becomes now independent of the random position of 
exfoliated flakes on a substrate and (2) the possibility to 
deposit different 2d materials on top of each other in order 
to build so-called van der Waals heterostructures (review 
about these heterostructures in reference ESI). 

It has been shown that these new transfer techniques 
can be used to overcome the problem of direct growth of 


electrode material onto graphene. For this the electrode 
structure is fabricated by means of e-beam lithography and 
metallization in a first step. Only afterwards the graphene 
flake is deposited on top of the prepatterned structure (fig¬ 
ure [lje)) |2[]. Therefore, the graphene flake is not exposed 
to the e-beam lithography step, hence resist residues and 
e-beam induced defects can be avoided. Furthermore, the 
growth conditions of the ferromagnetic electrodes and the 
oxide barrier are far more suitable and can be adapted to 
an adequate substrate. Especially, the whole expertise al¬ 
ready achieved in the fabrication of high quality magnetic 
tunnel junction devices (e.g. [TSHHU) may now also be ap¬ 
plied to graphene spin transport devices. In this respect, 
we especially refer to the improvement of magnetic tunnel 
junctions by annealing which results in crystallization of 
the interface between ferromagnetic metals and oxide bar¬ 
riers [80 [53]. The approach to first fabricate spin injection 
and detection electrodes on a wafer and then to deposit a 
stack of graphene on hexagonal boron nitride on top of it 
was e.g. applied in reference 0, where spin lifetimes of 
3.7 ns were measured at room temperature in trilayer gra¬ 
phene, which is the longest room temperature spin lifetime 
reported so far. Nevertheless, there might still be some is¬ 
sues with this new fabrication process because the insulat¬ 
ing barrier between graphene and ferromagnetic electrode 
is exposed to air at some point during the process. It is 
well known that many oxides [85] and particularly MgO 
[85H58] are hygroscopic. Therefore, it cannot be excluded 
that a partial hydroxylation of the MgO barrier may still 
limit the barrier quality in reference 0. 

As already mentioned, the possibility to deposit and 
stack various 2d materials on top of each other in order 
to fabricate van der Waals heterostructures is the second 
advantage of the transfer techniques. So far, many possi¬ 
ble 2d materials beyond graphene have emerged (overview 
e.g. in [ZSl EMH23) and combining these 2d materials with 
graphene has expanded the field of graphene research dra¬ 
matically. One of the first applications was the fabrication 
of stacks consisting of hexagonal boron nitride (hBN) and 
graphene. Due to the atomically smooth surface and sim¬ 
ilar lattice constant of hBN to graphene, these devices ex¬ 
hibit significantly improved electrical properties compared 
to graphene deposited on Si 02 [72] [931 IS]- The incorpo¬ 
ration of such hBN-graphene stacks in graphene spintronic 
devices led to a significant enhancement of the measured 
spin diffusion lengths as well as the charge carrier mobility 
Em. Another application of hBN which became feasible 
with the transfer techniques is its usage as an oxygen-free 
injection barrier when placed between graphene and the 
ferromagnetic electrodes [55] E2> Eg. By transferring exfo¬ 
liated hBN flakes, the above-mentioned problems in oxide 
barrier growth may be avoided. 

Finally, the transfer techniques offer a new route to 
functionalize graphene by putting different 2d materials 
in direct contact to graphene. One important applica¬ 
tion of such a heterostructure is the enhancement of the 
spin-orbit coupling (e.g. by the spin-orbit proximity effect 
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of tungsten disulfide on graphene 03) to allow for spin 
manipulation via the Bychkov-Rashba effect. We empha¬ 
size that many effects like the spin-orbit proximity effect 
greatly rely on contamination-free interfaces. However, de¬ 
pending on the exact kind of transfer technique different 
amounts of hydrocarbon contamination between the lay¬ 
ers of the heterostructures can be observed [951 [55] . Inter¬ 
estingly, there seems to be ” self-cleansing” effects in some 
heterostructures (e.g. graphene on hBN, M 0 S 2 , and WS 2 ), 
in which the surface contamination of the 2d materials au¬ 
tomatically aggregates into bubbles, leaving behind rather 
clean interfaces throughout ^m-sized graphene/hBN areas 

mm- 

3. Impact of metals on graphene 

In this section, we discuss the influence of metallic elec¬ 
trodes on both spin and charge properties of the under¬ 
lying graphene layer. One important tool to investigate 
the electronic properties is angle-resolved photoemission 
spectroscopy (ARPES), which can directly probe the band 
structure of graphene which is in direct contact to the 
metal (see e.g. references [501 EEHn!] for ARPES mea¬ 
surements of graphene in contact with Au, Ag, Fe, Ni, 
Co, Al, and Cu). Of special interest in the area of gra¬ 
phene spintronics are the 3d transition metals Fe, Ni, and 
Co, which are used as ferromagnetic electrodes to gener¬ 
ate spin polarized currents. ARPES measurements have 
revealed that all three metals significantly alter the band 
structure of graphene: The Dirac point of graphene shifts 
around 2.8 eV below the Fermi energy and the graphene’s 
7 r* band hybridizes with 3d bands of Co near the Fermi 

level mm- 

As a results, as soon as the ferromagnetic electrode gets 
into direct contact with graphene (which can be the case 
for pinholes in the separating oxide layer), the spin cur¬ 
rent is initially, i.e. right after electrical spin injection, no 
longer carried by pure graphene states near the Dirac point 
but rather by the hybridized states near the Fermi level. 
A significant higher spin scattering rate in such states in 
comparison to the states near the Dirac point of unmod¬ 
ified graphene may explain the overall short spin lifetime 
measured in devices with ohmic contacts exhibiting small 
contact resistance area products R C A < 1 kfl^m 2 [25] (see 
also figure [6]). Furthermore, the large density of states at 
the Fermi level due to hybridization results in a pinning of 
the Fermi level, which is equivalent to a screening of a gate 
electric field underneath the electrodes |1011fl03j . Consid¬ 
ering both the n-doping and the screening underneath the 
contacts caused by the Co, the gate voltage dependent 
doping profiles of figure [2jc) can be deduced for ohmic 
contacts, which in turn can explain the gate dependent 
resistance curve of such a device in figure [2]( a). 

A sufficiently thick oxide barrier without pinholes, which 
is preferable for spin transport devices, is expected to sup¬ 
press hybridization between the metal and graphene. Nev¬ 
ertheless, a field effect doping of graphene in such metal/- 


dielectric/graphene heterostructures is still expected be¬ 
cause of different respective work functions [10411106] . The 
thicker the insulating oxide barrier between graphene and 
metallic electrode (which results in larger R C A products) 
the lower the corresponding doping [106] , At some point 
the doping is small enough to be also tuned by the back 
gate voltage (see corresponding doping profile in figure 
[2jd) ). The different doping of graphene in between and un¬ 
derneath the electrodes can be seen by two distinct charge 
neutrality points (CNP) in gate dependent resistance mea¬ 
surements (figure [2jb)) : the right peak (blue arrow) corre¬ 
sponds to the graphene part between the electrodes (CNP 
mostly around a gate voltage of Vq = 0 V for untreated 
devices or at positive gate voltages because of p-doping 
during oxygen treatments, see reference E3) while the 
left one (red arrow) results from the graphene underneath 
the electrodes (CNP at large negative gate voltages be¬ 
cause of the n-doping of Co). For our Co/MgO/graphene 
devices, the contact-induced CNP at negative gate volt¬ 
age (see red arrow in figure [2](b) ) typically appears for 
R C A > lkfi/rm 2 [26] [27], which is also the value at which 
the I-V - curves of the contacts exhibit non-linear behav¬ 
ior m- At this point, we want to emphasize that only a 
sufficiently thick insulating barrier diminishes the metal- 
induced doping density of graphene strong enough that 
the contact-induced 2nd charge neutrality point becomes 
accessible by large back gate voltages. The required thick¬ 
ness depends on the used contact materials i 106 . There¬ 
fore, the value of R C A = 1 kfl/rm 2 is most likely only holds 
for our Co/MgO/graphene devices and may strongly differ 
for other material combinations. 

In the following, we discuss the origin of the contact- 
induced second CNP in more detail. Although this review 
article is about spin transport, we consider this discussion 
as crucial as it highlights two important aspects which are 
often neglected: (1) The significant modification of gra¬ 
phene’s transport properties right underneath the contacts 
and (2) the existence of lateral pn-junctions which the spin 
polarized charge carriers have to overcome. 

To confirm that the left CNP in our devices is indeed 
contact-induced, we now focus on a non-local spin-valve 
device with Co/MgO electrodes in more detail. This de¬ 
vice exhibit a large electrode-to-electrode variation of the 
respective R C A values. Figure [3][a) depicts the schematic 
layout. The measured contact resistance area products 
of the inner contacts 2 to 4 are: R C A 2 = 1.2kf2/.tm 2 , 
R c A 3 = 7.6kfl/im 2 , and R C A^ = 4.8kflp,m 2 . According 
to our previous studies [26] [27], we now assume that the 
interaction between Co and graphene is the dominating 
source of doping, whereas doping from fabrication-induced 
contaminations (like resist residues) or from the MgO it¬ 
self are subordinated. Following the above explanation 
that thicker barriers with larger R C A values diminishes Co¬ 
induced doping of graphene |104H106] . we can now assign 
larger doping levels to lower R C A values. This yields the 
potential profile of figure [3ja) for which we assume that 
the graphene between the contacts is hole doped by an ap- 
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Figure 3: (Color online) (a) Schematic layout of one device with assumed doping profile, which was deduced from the measured contact 
resistances (see text for more information), (b) and (d) show gate dependent resistance curves of regions A and B as a function of applied 
dc-current. (c) Gate voltage shift of the contact induced CNP normalized to the gate voltage position without applied dc-current as a function 
of dc-current. 


plied gate voltage. The depicted black curve corresponds 
to the position of graphene’s CNP relative to the Fermi 
level Ep (also see depicted Dirac cone). Gate dependent 
4-terminal resistance measurements were performed as a 
function of dc-current. For this, we apply a total current 
of I = / ac + 7dc = 1 /iA + I dc over the outermost electrodes 
1 and 5 and measure the voltage drop between contacts 2 
and 3 for region A and contacts 3 and 4 for region B by 
standard low frequency lock-in technique. 

As it can be seen for both regions in figures [3fb) and 
(d) both the position and magnitude of the left CNP sig¬ 
nificantly depends on the applied dc-current whereas for 
the right CNP of the bare graphene there is only a small 
decrease in resistance for both positive and negative dc 
current which can be explained by the larger local tem¬ 
peratures from current-induced Joule heating. Figure [3](c) 
depicts the dependence of the gate voltage position of the 
contact-induced left CNP as a function of dc-current. We 
show this data not only for regions A and B but also be¬ 
tween contacts 2 and 4, which we call region AB. Inter¬ 
estingly, the gate voltage shift of the CNP has different 
sign and amplitude in regions A and B. Both, the different 
sign and amplitude can also be seen in the doping pro¬ 
file of figure [3j a) (indicated by the arrows). Consistently, 
the dc-current dependent shift of the contact-induced CNP 
measured between contacts 2 and 4 (green curve in figure 


[3](c)) is in good agreement the arithmetic mean of the re¬ 
spective curves of regions A and B (blue curve in figure 

§c)). 

With the help of theoretical considerations about gra¬ 
phene field effect transistors |107tilll] we are able to qual¬ 
itatively understand the gate voltage shift of the contact- 
induced CNP. For this we assume that high quality tunnel 
barriers induce a held effect similar to top gates, which 
changes the Fermi level in graphene underneath the con¬ 
tacts. In fact, a reverse argumentation might be more intu¬ 
itive (we refer to reference Il06j ): Because of the different 
work functions between graphene and a metal a charge 
transfer will occur during the alignment of the Fermi lev¬ 
els. The charge transfer yields a carrier doping of gra¬ 
phene, but it also creates a voltage drop and therefore a 
local electric held over the separating dielectric MgO bar¬ 
rier as the barrier acts as a capacitor with the graphene 
and the metallic electrodes as its plates. These voltages 
Vtg* (TG as in top gate) are depicted in the schematic 
layout of a graphene held effect transistor in figure |4](a) . 
In the following, we assume that these voltages remain un¬ 
affected if the two inner contacts in figure |4](a) are used as 
voltage probes while the outer ones are source and drain 
contacts (this wiring applies to the measurement of hgure 

ID- 

Next, we determine the charge carrier density in the 
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Figure 4: (Color online)(a) Schematic layout of a graphene field effect transistor (see text for the definitions of the applied voltages). In case 
of the position marked with an X an equivalent circuit is depicted in (b), which can be used to calculate the charge carrier density at this 
point, (c) Illustration of the Fermi level shift by applying a bias voltage to both undoped and p-doped graphene, (d) Simplified explanation 
of the quantum capacitance of graphene (see text for more information), (e) Theoretical gate dependence of both the oxide capacitances and 
the graphene quantum capacitance. 


graphene flake. For the graphene underneath an inner 
electrode (marked by ”X” in figure |4](a)) we can apply the 
equivalent circuit in figure |4|(b). The MgO layer, which 
acts as a barrier for spin injection and spin detection, is 
approximated to be a top gate with capacitance Ctgii 
whereas the back gate voltage Vbg is applied over the ca¬ 
pacitance CbGi which represents the Si02 layer. Far more 
complicated is the incorporation of graphene in this pic¬ 
ture. For this we first discuss the Fermi level shift by 
applying a bias voltage over the graphene flake. In figure 
Qc) this effect is illustrated in a very simplified system, in 
which e.g. the doping effects of the contacts are neglected. 
We first consider ideal, undoped graphene, in which case 
the Fermi level of the graphene underneath the grounded 
contact (left contact in figure |4](c)) is at the Dirac point. 
A bias voltage applied to the right contact results in a 
lateral voltage gradient along the graphene channel with 
corresponding shifts of the Fermi level fwiirns] , Along 
the whole graphene flake a gradual nn’-junction is gener¬ 
ated, which can be verified by photocurrent measurements 

m- 

The Fermi level shift underneath the contacts with 
the applied source-drain-voltage Vsd can be written as 
Ep = — eV c h where the energy scale is set to zero at the 
CNP of graphene. In this simplified system the voltage 
drop V c h along the graphene channel only depends on the 
electrostatic potential by the applied bias. But in gen¬ 
eral, the Fermi level also depends on the electrochemi¬ 


cal potential. Accordingly, the position of the Fermi level 
Ep = —eV c h also depends on other factors such as the 
applied gate voltage or changes in the electrochemical po¬ 
tential by adsorbates. This is also illustrated in figure Qc) 
in the case that graphene is p-doped by an applied gate 
voltage or adsorbates. Apparently, the local area with 
vanishing charge carrier density, i.e. where the Fermi level 
is at the Dirac point, moves along the graphene channel 
(from left to right in the lower panel of figureQc)) with in¬ 
creasing p-doping concentration. Such a shift of the CNP 
by an applied gate voltage is e.g. observed in photocurrent 
and thermal infrared microscopy [112H115] . 

The capacitance of the graphene part in the equivalent 
circuit over which the voltage Wh drops is the quantum 
capacitance Cq of graphene 1l’fi 117] . The quantum ca¬ 
pacitance directly results from the Pauli principle and be¬ 
comes relevant in materials with a small density of states 
(DOS) near the Fermi level. To explain the quantum ca¬ 
pacitance, we first consider ideal graphene with the Fermi 
level at the Dirac point, where the DOS vanishes (figure 
B d )> left band structure). If now the charge carrier den¬ 
sity in graphene is changed by An the Fermi level shifts 
by AEpi as the electrons have to occupy higher energy 
states. For single layer graphene, the Fermi level Ep and 
the DOS at the Fermi level D(Ep) are given by [1181 : 

Ep = hvpy/Trn , D(Ep) = 2 ^ F 2 , (1) 

7r (Tiv f) 
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where v F is the Fermi velocity. If now the charge car¬ 
rier density is doubled by adding the same amount of 
charges An a second time, the increase of the Fermi level 
becomes less strong because of the previous increase in 
DOS (AEp 2 < 2Ai?Fi) (see figure Qd)). As mentioned 
above, the electrons at the Fermi level have the energy 
Ey = —eV c h . Combining this expression with equation 0 
and the general definition of a capacitance Cq = dQ/dV, = h 
(with the charge Q = Ne = Ane and area A) results in 
the quantum capacitance of graphene: 

A 7T (hv F ) n \hv F J y/irhv F 

In figure |4](e) we plot this quantum capacitance of gra¬ 
phene as a function of the applied gate voltage (black solid 
line). Due to broadening of the DOS by both thermal ex¬ 
citation at room temperature and by structural inhomo¬ 
geneities the experimentally determined graphene quan¬ 
tum capacitance is also broadened [119( H2D] which is de¬ 
picted as a black dashed line. Furthermore, we include 
the capacitance of both oxide layers (the MgO of the top 
electrode with thicknesses of 2 and 3 nm and the Si02 of 
the back gate). 

Now we go back to the equivalent circuit in figure[4](b). 
First, we start with the graphene part between the con¬ 
tacts (marked by a circle in figure |4](a)). As there is no top 
gate capacitance in this region, back gate and quantum ca¬ 
pacitances are put in series and give a total capacitance 
Ctotai = CbgC'q/(C'bg + Cq)- As the back gate capaci¬ 
tance is much smaller than the quantum capacitance (com¬ 
pare to figure (4^e)) , it dominates the total capacitance: 
Ctotai ~ Cbg- Accordingly, the quantum capacitance does 
not play any significant role in the back gate induced field 
effect (see also )116l 1117] ). This is an important conclu¬ 
sion as the bias voltage can only change the charge carrier 
density at the node of the equivalent circuit in figure |4])b) 
with help of the quantum capacitance. The fact that the 
quantum capacitance can be neglected now explains why 
the CNP of the bare graphene part between the contacts 
does not shift with dc currents (figures [3](b) and (d)). In 
contrast, the quantum capacitance plays a significant role 
for graphene parts underneath the electrodes as it is in 
the same order of magnitude as the electrode capacitance 
(figure (4je)). If we assume a fixed charge carrier density 
at the node of the equivalent circuit of figure |4])b) for the 
contact region, now the shift of the contact-induced CNP 
in figures [3])b) and (d) becomes clear: By changing the dc 
current (or accordingly the bias voltage) the voltage V c h 
will change according to the discussion in figure [d|c). If 
we want to keep the charge carrier density at the node of 
the equivalent circuit at n = 0 (i.e. to the contact-induced 
CNP) the voltages Vtgi (be. the actual doping of the gra¬ 
phene by the contacts) and Vbg also have to change. To 
simulate this, all changes of Vtgii Vbg> and Wh have to 
be calculated self-consistently (e.g. a change in the back 
gate voltage will shift the Fermi level and therefore the 


value of the graphene quantum capacitance), which goes 
beyond the scope of this paper. For a far more detailed 
description we refer to the references DEZHm]. 

Finally, we again point to the fact that the different 
charge transport properties of graphene underneath the 
contacts and graphene in between the contacts are so far 
widely neglected. It has to be seen how the change in the 
doping underneath the contacts with different dc-biases 
or the pn-junctions along the graphene channel have an 
impact on the performance of graphene-based spintronic 
devices fT2TlH2^| . 

4. Determination of spin lifetimes 

In this section we discuss on how to extract the spin 
lifetime out of spin precession measurements and address 
how several limitations to this approach have recently been 
recognized. We restrict ourselves to non-local spin trans¬ 
port measurements, so we can neglect the drift term in 
the Bloch-Torrey equation and only consider changes in 
the net spin vector s by spin precession about a perpen¬ 
dicular magnetic field B (see measurement configuration 
in figure a)), spin diffusion, and spin dephasing and re¬ 
laxation [2]: 

g = sx^ 0 + D s V 2 s~- . (3) 

at t s 

Here cJq = g^sB/U is the Larmor frequency, g is the gyro- 
magnetic ratio, D s the spin diffusion constant, and r s the 
spin lifetime. 

The advantage of non-local measurements is that the 
spin signal can be significantly decoupled from spurious 
charge signals [2 [3l [123] . To achieve such a non-local 
configuration, dc measurements or lock-in techniques with 
low frequencies are needed, because electric pulses in case 
of RF measurements propagate through the whole device, 
which means that no non-local part can exist. The dc or 
low frequency measurements lead to a stationary or quasi¬ 
stationary condition and therefore the time derivative of 
equation © is set to zero. Hence a solution to 

sxi3o + AV 2 s-- =0. (4) 

Ts 

has to be found to extract the spin lifetime by means of a 
fit to the dc spin precession curve (so-called Hanle curve). 
Typical Hanle curves for a spin transport device are shown 
in figure [5](b) with both parallel and antiparallel align¬ 
ments of the inner Co electrodes in figure [5ja) . The non¬ 
local spin resistance AI? n i can be determined at B = 0. 
As expected, both Hanle curves merge at larger magnetic 
fields. However, they do not become constant but rather 
increase above \B\ > 0.2 T. We note that this magnetic 
field dependent background signal is typical for most stud¬ 
ies and it can even be much more pronounced. One ex¬ 
planation for this background is that the magnetization of 
the ferromagnetic electrodes can rotate out-of-plane with 
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Figure 5: (Color online) (a) Schematic cross section with wiring con¬ 
figuration for non-local spin-transport measurements, (b) Hanle spin 
precession measurement of a bilayer graphene device for a perpendic¬ 
ular magnetic field sweep with parallel and antiparallel alignments 
of the respective spin injection and detection electrodes 

increasing perpendicular magnetic field [0} 112414126] . In 
principle, this rotation of the magnetization can account 
for a background signal which is symmetric in magnetic 
field. But we often observe an antisymmetric background 
signal with a linear term in B , which also dependents on 
the charge carrier density and the wiring of the device 
m- These findings cannot be explained by the rotation 
of the electrode’s magnetization alone. We have recently 
shown that such a background of second polynomial order 
can be caused by an inhomogeneous current flow through 
the oxide barriers. As a result there is a charge accumula¬ 
tion signal next to the actual spin accumulation signal in 
the non-local voltage which can be explained by a redis¬ 
tribution of charge carriers by a perpendicular magnetic 
field similar to the classical Hall effect m- 

One fundamental drawback of every solution to equa¬ 
tion 0 is that the three parameters g, D s , and r s can¬ 
not be determined independently. This is because of the 
steady state condition, which allows the multiplication of 
equation Q by any factor a without changing the over¬ 
all result. Therefore a fit using a solution to equation @ 
is invariant with respect to the transition (g,D s ,r s ) —» 
( ag , qD s , T s /a). Hence one of the three parameters has to 
be assumed or determined by other measurements. 

Often g = 2 is assumed for the analysis of spin preces¬ 
sion measurements. For pristine or moderately modified 


graphene flakes this assumption is confirmed by electron 
spin resonance (ESR) measurements 02] 03} ll2S] . It re¬ 
mains to be seen if this value also holds for functional¬ 
ized graphene. By now, spin precession measurements of 
hydrogenated graphene suggest larger ^-factors [12914131] . 
Of course, the most elegant way to circumvent this uncer¬ 
tainty is to give up on the non-local, steady state condition 
and apply time-resolved RF measurements which offer the 
possibility to directly determine the g-factor from time de¬ 
pendent spin precession. 

In the following, we discuss some assumptions and sim¬ 
plifications, which have to be made in order to find useable 
solutions to equation Q. One of the first and still widely 
used analytical solution was developed by Johnson and 
Silsbee 023]. But in recent years it was found that the 
measured spin lifetime of a graphene device scales with 
the R C A product at least up to values of several tens of 
kRLgm 2 (see figure [0) [20} [20} [27} 03] ) suggesting that the 
measured spin lifetime is not the intrinsic spin lifetime of 
graphene but rather has extrinsic origin. One route to ex¬ 
plain this dependence is that the measured spin lifetime 
can be significantly underestimated if spin relaxation by 
the contacts is not included in the fit model and hence 
some work was done in the direction of more elaborated 
models [55 ]fT321fT33] . But even these models make two 
crucial assumptions: (1) The injection barriers are homo¬ 
geneous and can be characterized by the R C A value only 
and (2) the spin lifetime in graphene underneath the con¬ 
tacts is the same as the spin lifetime in the bare graphene 
part between the contacts. 

As it was already discussed in section [3] the growth of 
homogeneous oxide barriers on graphene is quite challeng¬ 
ing. But as long as the contact is not spatially homoge¬ 
nous, the measured contact resistance can only be an aver¬ 
aged value. Especially if there are pinholes within the bar¬ 
rier the contact resistance area product cannot precisely 
be determined, because the exact number and sizes of the 
pinholes cannot be deduced from electrical measurements. 
Therefore, the uncertainty in the contact resistance will 
also lead to an uncertainty in the extracted spin lifetimes 
by the aforementioned models. 

And also the second assumption of a single spin lifetime 
seems to be oversimplified considering section [2] As it was 
demonstrated by ARPES measurements 11 (')(') the in¬ 
teraction between graphene and ferromagnetic metals can 
significantly change the band structure of graphene due to 
hybridization. The hybridized states are directly at the 
Fermi level and therefore both charge and spin transport 
will occur through these states, which may have a strong 
impact on the spin lifetime. But even if we assume a suf¬ 
ficiently thick oxide barrier, which inhibits this hybridiza¬ 
tion, the assumption of a single spin lifetime may still be 
too simple if the spin lifetime depends on the Fermi level 
or, correspondingly, the charge carrier density. Because 
the field effect doping of the contacts }l06j can e.g. lead to 
the case at which the graphene underneath the contacts is 
highly n-doped whereas the graphene between the contacts 
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Figure 6: (Color online) Spin lifetime vs averaged contact resistance 
area product R C A of respective injection and detection electrode at 
an electron density of n = 1.5 X 10 12 1/cm 2 at room temperature. 
Shown are data from our old fabrication method (as-fabricated single 
and bilayer graphene devices from reference m) and three devices 
where the contact resistances were successive increased by oxygen 
treatments from reference [27]) as well as data from our new bottom- 
up fabrication method from reference (Bj. 

is at its CNP. By now, only some research was done on the 
topic how two spatially different graphene parts with dif¬ 
ferent spin relaxation times influence the overall measured 
spin lifetime (e.g. in reference [15|). 

Finally, we note that a rough MgO layer like the one 
in figure [lid) not only can cause pinholes but also leads to 
an equally rough surface of the Co electrode which is de¬ 
posited on top. But such a rough ferromagnetic interface 
can yield stray fields, which may result in a broadening of 
the spin precession curve and thereby an apparent reduc¬ 
tion of the extracted spin lifetime [134 . 

5. Determination of charge carrier mobilities 

In this section we discuss the determination of the 
charge carrier mobility /i from gate dependent resistance 
measurements. There are two reasons making this an im¬ 
portant topic for graphene-based spin transport studies. 
Firstly, the dependence of the spin lifetime and the charge 
carrier mobility is often used to identify the dominating 
spin relaxation mechanism mm- Secondly, the use of ox¬ 
ide barriers for spin injection and detection can lead to 
the appearance of the contact-induced CNP as it has been 
discussed in section [3] We will demonstrate that this 2nd 
CNP has significant impact on the determination of the 
charge carrier mobility. 

Firstly, we analyze the electron mobility in one of our 
devices by different methods which are currently used in 
literature showing that the determination of the mobility 


is not unambiguous and yield values which vary by more 
than a factor of 2. The simplest model considers a gra¬ 
phene field effect transistor, for which a linear dependence 
between the charge carrier density n (or accordingly back 
gate voltage) and the conductivity a is often observed (see 
e.g. [l35j ). For this case several groups (e.g. [T3BHISSI) as¬ 
sume a simple Drude model and define the Drude mobility 
as: 

MDrude = a-/(ne). (5) 

We emphasize that this Drude mobility is extracted at 
a single point of the conductivity curve. To account for 
the residual conductivity a o even in the case of n —» 0 the 
conductivity can be written as a = )ien + <7o- Starting 
from this equation the extracted field effect mobility /ife 
from fitting the slope of the conductivity curves is given 
by (e.g. in references 52, 55]): 

/^fe = (1 /e) (Act/A n). (6) 

We applied both methods to a conductivity curve from 
one of our actual spin valve devices in figure [7](b) . The red 
circle marks the point at which le and Mfe have been 
determined. Remarkably, Drude exceeds /r fe by almost 
60%. In the as-fabricated state of this device the conduc¬ 
tivity drops only slightly at large negative gate voltages, 
which indicates that the contact induced CNP is close but 
nevertheless out of range of the applicable gate voltage. 
After an oxygen treatment which primarily reduces the 
doping effect of the electrodes m the contact-induced 
CNP shifts into the accessible gate voltage range (figure 
[?Je)) and furthermore complicates the analysis of the car¬ 
rier mobilities. 

This can be seen in the resistance curve of another 
spin valve device in figure [7)(a) that also shows two distinct 
CNPs. We can decompose resistance contributions of both 
graphene regions (underneath and in between the con¬ 
tacts; model for fitting is discussed further below), which 
is shown by red and blue dashed lines, respectively. The 
overlapping of both curves clearly demonstrates one major 
issue of the two mobility fits discussed so far: As soon as 
the doping underneath the contact approaches the one of 
the bare graphene part, both the absolute value and the 
slope of the measured gate dependent resistance curve sig¬ 
nificantly changes in the overlapping part. And so does the 
extracted mobility. The reason for this is quite clear: Both 
equation (|5j) and ([6]) are based on one well-defined trans¬ 
port regime, whereas the contact-induced CNP points to 
the fact that the charge transport in graphene underneath 
and in between the contacts can differ significantly. 

But before we discuss more detailed fitting models, 
which include different transport regimes, we first want 
to discuss another issue when using equation (Jiij) which 
arises as soon as the conductivity curves are not linear as 
the ones in figure [7j In this case, the extracted mobility 
/i, which is determined by the slope, depend on the charge 
carrier density n. But this is a mathematical contradiction 
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Figure 7: (Color online) (a) Decomposition of a gate dependent resistance measurement into two Dirac curves for device A corresponding 
to the n-doped graphene underneath the contacts and the undoped graphene between the contacts. Figures (b)-(d) show the conductivity 
curve of device B, which is fitted with different models for the determination of the charge carrier mobility (dashed lines). See text for more 
information, (e) After the contacts of the same device B were oxidized, the contact induced CNP shifted into the measurable gate voltage 
range. 


which is often neglected. The derivative of a with respect 
to n is now given by: 


dcr 

dn 

=> pip) 


— (nep(n) + a 0 ) 
an 
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break down near the CNP. But there are also models try¬ 
ing to fit the whole conductivity curve including the two 
CNPs (e.g. in reference HU). We also use the following 
model which assumes two differently doped graphene parts 
with different mobilities /xp 


1 

a 


E 


i 

V O? + nf 0 ) e/xi 


+ Po, 


( 10 ) 


The term d/x(n)/dn is neglected by the approach of 
equation ([6]) and only vanishes for strictly linear conduc¬ 
tivity curves. 

But the question arises if the mobility indeed depends 
on the charge carrier density in case of a non-linear con¬ 
ductivity curve as the non-linearity can also be explained 
by a constant mobility (d/x/dn = 0) if there is a charge car¬ 
rier independent contribution po to the overall resistivity. 
Such a contribution can result from both short- and corre¬ 
lated long-range disorder [139 ;. Hence, the corresponding 
conductivity is given by (see e.g. references [7H 1'40]) : 

1 1 

- — -—-1 Po- (9) 

a nep * + a o 

Figure [tJc) demonstrates that a significant part of the 
electron branch for Vq > 0 of the non-linear conductivity 
curve can be fitted under the assumption of the constant 
mobility p* = 8380cm 2 /Vs. However, the fitting fails com¬ 
pletely in the hole branch (not shown), because here the 
impact of the contact-induced CNP has to be considered. 
Furthermore, all three approaches that are discussed so far 


where n, = a(Vc, — Vcnp.i) is the gate voltage Vq de¬ 
pendent charge carrier density, a = 7.18 x lO 10 l/(cm 2 V) 
the capacitive coupling for a 300 nm thick SiC >2 layer, n*o 
the residual charge carrier densities due to the presence of 
electron-hole puddles [142j and thermally excited carriers, 
which prevent a divergence at the CNPs [ 141] , and a gate 
independent resistivity po. This model can be fitted quite 
well to the conductivity curve of figure [7](d) (the given mo¬ 
bility p = 10900 cnr 2 /Vs corresponds to the Dirac curve 
with its CNP near zero gate voltage). 

But now there are two issues: At first, the values of the 
charge carrier mobilities as extracted from the four mod¬ 
els spread over a broad range between 4580cm 2 /(Vs) and 
10900 cm 2 /(Vs). Secondly, there is a serious issue when 
extracting the mobility from equation (10), which already 
contains seven parameters, as it cannot be used to ade¬ 
quately fit the conductivity curve of the same device af¬ 
ter oxygen treatment (figure [7](e)). We tried several other 
models to fit the curve in figure [7](e) , but for a reasonable 
fit we need at least eight free parameters. But as we al¬ 
ready briefly discussed in one of our previous publications 
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[261 and will discuss in this section in more detail, such a 
large number of free parameters make such a model use¬ 
less. Because as long as the Dirac curves have symmetric 
shapes, which they have in figure [7ja) and in our previous 
publications [261 E7], only seven parameters can reliably 
be extracted out of a gate dependent resistivity curve with 
two visible CNPs: The position, magnitude and width for 
each of the two peaks and the background. Therefore, a 
model of eight parameters is overdetermined and accord¬ 
ingly we are able to change the mobility values in such a 
model over an unreasonable wide range without seriously 
worsening the fit results. 

To approach reliable mobility fits, we now deduce sev¬ 
eral parameters which have to be included into a conduc¬ 
tivity model to adequately fit a gate dependent resistivity 
curve with two CNPs. First, there is the carrier mobility 
as the actual quantity of interest. As already discussed in 
section [3] the interaction with the electrodes can signifi¬ 
cantly change the electronic properties of graphene. We 
therefore have to consider two mobilities for each device: 
one for the graphene underneath the contacts and one for 
the graphene area in between the contacts. The same ar¬ 
gument holds for the electrochemical doping of the respec¬ 
tive graphene parts, which is obvious as we observe two 
CNPs. Two additional parameters are needed to include 
spatial variations of the electrochemical potentials within 
each graphene region. These two parameters can also ac¬ 
count for the minimum conductivities as the variations 
in the potential lead to the electron-hole puddles at the 
charge neutrality point I l4ll . So far this yields six param¬ 
eters. 

In the next step, we have to critically review the gate 
voltage as the tuning parameter for the resistivity mea¬ 
surement. As discussed e.g. in | il01H103j the interaction 
between contact and graphene may lead to a screening of 
the gate electric field underneath the electrodes. We as¬ 
sume that there is a gradual transition from pinning to 
depinning of the Fermi level underneath the contacts with 
increasing oxide barrier thickness and quality. Therefore, 
we have to consider a screening factor as a seventh pa¬ 
rameter which accounts for the effective Fermi level shift 
in the graphene underneath the contacts as a function of 
applied gate voltage. 

Furthermore, there is the unknown transition of the 
electrochemical potentials between contact covered and 
bare graphene part, i.e. the exact shape and lateral ex¬ 
tension of the pn-junctions, which are known to exist near 
the edges of the contacts [102 . 103 L In the most sim¬ 
ple approximation, at least one more parameter has to be 
considered which describes the decaying length of the pn- 
junction (see e.g. references (1431 1144] b And finally, there 
might be a gate voltage independent contribution po to 
the overall resistance as already discussed in equation Q. 
Hence, at least nine parameters are necessary to simulate 
the gate dependent resistivity, whereas only seven indepen¬ 
dent parameters may be experimentally determined by a 
gate dependent resistance measurement. 


The whole problem with the analysis of the carrier mo¬ 
bility from gate dependent resistance measurements is that 
both the contact-covered and bare graphene parts are mea¬ 
sured simultaneously in series. The lack of independently 
probing both regions is bound to the very nature of a pure 
electrical transport measurement. Therefore, independent 
measurements have to be conducted to determine some of 
the aforementioned parameters separately. Especially, in¬ 
formation about the spatial change in the electrochemical 
potential along the contact induced pn-junctions may be 
helpful. These can be obtained from scanning photocur¬ 
rent microscopy [1021110311145] . 

Finally, we note that most other graphene-based nano- 
electronic devices are fabricated with low-ohnric, metallic 
contacts. This might be the reason why the aforemen¬ 
tioned issues concerning the determination of the charge 
carrier mobility is not discussed in recent publications in 
more detail. For metallic contacts, a complete pinning of 
the Fermi level underneath the contacts can be assumed 
(see section [3]). Therefore, the complete contact area of 
the device only contributes as a gate independent con¬ 
stant resistance to the overall gate dependent resistance 
(e.g. p 0 in equation ([ 9 ])). Additionally, the impact of both 
the contact area and the pn-junctions on the transport 
measurement decreases with increasing separation of the 
electrodes. We note that in many charge-based electri¬ 
cal transport studies the contact separation is much larger 
than in graphene-based spin transport devices. 

6. Critical review and perspectives 

In this final section, we critically review some exper¬ 
imental studies which were used to investigate the domi¬ 
nant spin relaxation mechanism in graphene. Maybe the 
most direct way to identify the dominating spin relax¬ 
ation mechanism is to evaluate the dependence of the spin 
lifetime on the charge carrier mobility 0E[. A linear 
dependence of r s on p or t p , which is the momentum 
scattering time, is a priori suggestive of an Elliott-Yafet 
(EY) spin scattering mechanism, while the inverse relation 
t s oc 1/p oc 1 /t p indicate the dominance of a D’yakonov- 
Perel’-like (DP) spin scattering mechanism. In our previ¬ 
ous studies we found an inverse dependence of spin lifetime 
on the mobility for both single-layer [26| and bilayer [2S] 
graphene which is thus indicative for DP-like spin scatter¬ 
ing. This dependence was only found in devices with con¬ 
tact resistance area products larger than R C A > 1 kflpm 2 
(figure pi a)). In contrast, all devices with R C A < 1 kttpm 2 
(figure 481(a)) and devices where the contact resistances 
were enhanced by subsequent oxygen treatments (figure 
[8](b)) do not show this l//r-dependence. All these devices 
have been prepared by the conventional top-down method 
for which the MgO injection and detection barrier was di¬ 
rectly deposited onto graphene. 

Most strikingly, the new generation of devices where 
we apply the bottom-up approach by transferring a hBN- 
graphene-stack on prepatterned electrodes (see figure [TJe) 
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Figure 8: (Color online) (a) Spin lifetime vs electron mobility at n = 1.5 X 10 12 1/cm 2 at room temperature for as-fabricated single and 
bilayer graphene devices. Only for devices with R C A > 1 kfl//m 2 a D’yakonov-Perel’ like l//x-dependence can be fitted to the data. The 
line for the devices with R C A < 1 kQ/rm 2 is just a guide to the eye and illustrates the deviation from a 1///-dependence, (b) Same trends 
as in figure (a) combined with the development of the spin lifetime and mobility of four devices during oxygen treatment, (c) Devices with 
R C A > 1 kfl/zm 2 both fabricated with our old process (exfoliated graphene on Si02 with subsequent deposition of the electrodes) and the 
new one (prepatterned electrodes and transfer of a hBN/graphene-stack like depicted in figure Qe)). Data taken from references [5l 1261 [271 . 


and section [2| exhibit significantly enhanced charge and 
spin transport properties (see full symbols in figure |8|[c)). 
We attribute the increase in mobility to the hBN substrate 
while we relate the increase in spin lifetime to improved 
contact properties according to our advanced transfer tech¬ 
nique which has several advantages over the previous fab¬ 
rication methods. At first, the contact region has not been 
exposed to an electron beam which may reduce the num¬ 
ber of defects in graphene. [57] [58] Secondly, the interface 
between MgO and graphene is expected to be of higher 
quality yielding more homogeneous barriers which can be 
seen by the larger R C A values for devices with thinner 
MgO layer thicknesses [5]. 

Our results indicate that the overall improvement of 
the spin properties primarily result from the improvement 
of the contact properties suggesting that the observed 1 //x 
dependence in the initial work is of extrinsic origin. In this 
context, we again want to emphasize (sections [4] and [5]) 
that the determination of both the charge carrier mobility 
and the spin lifetime is by no means unambiguous and may 
thus result in contradicting r s — t p dependencies. 

Furthermore, we demonstrated in (26] and m that 
the spin lifetimes in |25j are most likely limited by con¬ 
tact properties (also see figure [6| in all devices. In fact, 
there are indications that even in our newest bottom-up 
devices the contacts may still limit spin transport proper¬ 
ties (paper in preparation). If in current devices extrinsic 
parameters limit spin transport properties, it will be inter¬ 
esting to see how far more advanced fabrication methods 
will yield devices with even longer spin lifetimes and larger 
carrier mobilities which ultimately allows to unveil intrin¬ 
sic spin scattering mechanisms. 
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